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We show that the lightest right-handed sneutrino in TeV scale supersymmetric B — L model with 
inverse seesaw mechanism is a viable candidate for cold dark matter. We find that it accounts for 
the observed dark matter relic abundance in a wide range of parameter space. The spin-independent 
cross section oi B — L right-handed sneutrino is consistent with the recent results CDMS II and 
XENON experiments and it is detectable in future direct detection experiments. Although the B — L 
right-handed sneutrinos annihilate into leptons, the PAMELA results can not be explained in this 
model unless a huge boost factor is considered. Also the muon flux generated by B — L right-handed 
sneutrino in the galactic center is smaller than Super-Kamiokande's upper bound. 

PACS numbers: 



I. INTRODUCTION 

The experimental verifications of non-vanishing neutrino masses and the alluring hints of dark matter's 
(DM's) existence are serious indications for new physics beyond the Standard Model (SM). Supersymmetry 
(SUSY) is an attractive candidate for new physics at TeV scale that provides an elegent solution for the 
SM gauge hierarchy problem and stabilize the SM Higgs mass at the electroweak scale. The minimal 
supersymmetric standard model (MSSM) is the simplest extension of the SM. In order to account for the 
observed neutrino masses and mixing, SM singlets (right-handed neutrinos) are usually introduced. 

In the MSSM with i?-parity conservation and universal soft SUSY breaking terms, the lightest neutralino, 
which is typically bino dominated, is an attractive candidate for cold DM. However, the current experimental 
constraints on SUSY particles lead to overproduction of bino relic abundance, in contradiction with the 
observational limits of the Wilkinson Microwave Anisotropy Probe (WMAP) In addition, the recent 
results of Cryogenic DM Search (CDMS II) [S] set an upper limit on the DM-nucleon elastic scattering spin 
independent cross section of order 3.8 x 10~''^cm^ for DM mass of 70 GeV, which imposes stringent limits 
on the lightest neutralino of the MSSM, even if it consists of gaugino-Higgsino mixture. Therefore, one 
concludes that the DM constraints severely reduce the allowed range in the parameter space of the MSSM. 
It is also worth mentioning that the observed anomalies in the cosmic rays may favor the type of dark matter 
that annihilates into leptons not to quarks, unlike the lightest neutrlino in MSSM. 

TeV scale right-handed neutrinos can be naturally implemented in supersymmeric B — L extension of the 
SM (SUSY B - L), which is based on the gauge group Gb-l = 5'J7(3)c x SU{2)l x U{1)y x U{1)b-l 
In this model, three SM singlet fermions arise quite naturally due to the U{1)b-l anomaly cancellation 
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conditions. These particles are accounted for right-handed neutrinos, and hence a natural explanation for 
the seesaw mechanism is obtained 0-0]. This work was based on earlier papers in Ref. 01. In this class of 
model, the scale of B — L symmetry breaking is related to supersymmetry breaking scale [7[. Therefore, the 
right-handed neutrino masses are naturally of order TeV scale. This has initiated a considerable interest in 
analyzing the phenomenological implications of these models and their possible signatures at the LHC Q . 

TeV scale SUSY B — L \s one of the simplest models that provides viable and testable solution to the two 
puzzles of the DM and the neutrino masses. However, in order to fulfill the experimental measurements for 
the light-neutrino masses, one of the following scenarios must be adopted: {i) Type I seesaw mechanism with 
very small Dirac neutrino Yukawa couplings, Yi, < 0(10^^) {ii) Inverse seesaw mechanism with order 
one Yukawa couplings and small mass scale ^ C>{1) KcV, that corresponds to scale of breaking a remnant 
discrete symmetry, (—1)^"'"'^ This work was based on earlier papers in Ref. [11]. In the first case, due to 
the smallness of Dirac Yukawa couplings, the right-handed neutrino sector has a very suppressed interaction 
with the SM particle. Therefore, the prediction of SUSY B — L remains close to the MSSM ones. It turns 



out that the DM candidate of this model is still the lightest neutralino 1^ , which is a kind of a mixture of 
three neutral gauginos and four neutral Higgsino. 

In this paper we consider the scenario where the right-handed sneutrino in SUSY _B — L with inverse 
seesaw is the lightest SUSY particle (LSP) and stable, so that it can be a cold DM candidate lia]. It is 
worth mentioning that in MSSM the left-handed sneutrino is the only stable weakly coupled neutral boson 
that can be a DM candidate IJ]. However, the current limits in direct detection experiments rule out 
this possibility, since left-handed sneutrino has a tree level interaction with the Z gauge boson, hence its 
elastic scattering cross section with nucleons is quite large. In case of MSSM extension with TeV right-handed 
neutrino (N) superfields, the interaction of N with the SM particles can be obtained from the superpotential: 

W^Wmssm + Y,N'LH2. (1) 

Due to smallness of YJy, the annihilation cross section of the right-handed sneutrinos is very suppressed, 
hence its relic density is larger than the WMAP limit. Moreover, since they cannot couple to the quark 
sectors, the direct detection experiments such as CDMS II cannot be tested in this case. 

We will show that in the SUSY B — L model with Y^, ^ 0{1) the dominant annihilation channel for right- 
handed sneutrino is given by the four point interactions, leading to and /i". While the effective couplings 
of right-handed sneutrinos with quarks are obtained through the exchange of TeV scale B — L gauge boson, 
Zb^l- Thus, one finds that the WMAP results of the DM relic abundance and CDMS II/XENON for direct 
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detection can be accommodated in SUSY B — L simultaneously 

The paper is organized as follows. In section 2 we analyze the supersymmetric B — L model with Inverse 
Seesaw Mechanism. We study the neutrino sector and the neutrino mass eigenstates. In section 3 we analyze 
the B — L right-handed sneutrino mass and interactions. Section 4 is devoted for the analysis of the relic 
abundance oi B — L right-handed sneutrino. We show that the WMAP limits can be easily satisfied in a wide 
range of parameter space. In section 5 we discuss the direct detection rate oi B — L right-sneutrino DM. We 
show that the elastic cross section of our DM candidate with nucleon is consistent with the recent results of 
CDMS II/XENON experiment and it is detectable in near future experiments. In section 6 we discuss the 
indirect detection rate oi B — L right-handed sneutrino. We show that the annihilation channels of right- 
handed sneutrino into leptons are subdominant, therefore it cannot account for the controversial PAMELA 
results. We also show that the muon flux generated from the right-handed sneutrino in the galactic center 
is much smaller than the Super-Kamiokande's limits. Finally we give our conclusions in section 6. 
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II. SUPERSYMMETRIC L MODEL WITH INVERSE SEESAW MECHANISM 

The proposed TeV scale supersymmetric B — L extension of the SM is based on the gauge group SU (3)c x 
SU{2)l X U{1)y X U{1)b-l, where the U{1)b-l is spontaneously broken by a chiral singlet superfield xi 
with B~ L charge = +1 and X2 with B — L charge = — 1. As in the conventional B — L model, a gauge boson 
Zb-l and three chiral singlet sueperfields Ni with B — L charge = — 1 arc introduced for the consistency 
of the model. Finally, three chiral singlet superfields Si with B — L charge = +2 and three chiral singlet 
superfields S2 with B — L charge = — 2 are considered to implement the inverse seesaw mechanism. The 
superpotential of the leptonic sector of this model is given by 

W - Y.E'LHi + Y,N'LH2 + YsN^xiS2 + ^^^1^2 + m'xiX2. (2) 

It is worth noting that the chiral singlet superfields X2 and N have the same B — L charge. Therefore, one 
may impose a discrete symmetry in order to distinguish them and to prohibit other terms beyond those given 
in Eq. In this case, the relevant soft SUSY breaking terms, assuming the usual universality assumptions, 
are as follows 

- Csoft = "^^I-^I' + YJ^N'LH2 + Yj^E'LHi + Y/iV^^aXi + B^iHiH2 + B/i'xiX2 

+ ^MiBB +h'l2W''W'' + ^A'hrr + lMB-LZB-LZB-L + h.c, (3) 



where the sum in the first term runs over (j) = Hi,H2,Xi,X2,L, E'^, N'^, Si, S2 and Y^ = YlAl {L = e, S) 
is the trilinear associated with lepton Yukawa coupling. In order to prohibit a possible large mass term 
MS1S2 in the above, we assume that the particles, N^, Xi,2, and S2 are even under matter parity, while 
is an odd particle. The B — L symmetry is radiatively broken by the non-vanishing vacuume expectation 
values (VEVs) (xi) — v[ and (X2) = ^2 H- The tree level potential F(xi,X2) is given by 

^(X1,X2) = ^5l-L(IX2p - IxiH' + A^^IXiP + M2IX2I' - /i3(XiX2 + h.c). (4) 

At GUT scale, fif = TOq + fi'^ ,i — 1,2 and jJ^ = -Bji'. However, they have different evolution from GUT 
scale to TeV scale and ^2 becomes negative, so that B — L is spontaneously broken fvl]. The minimization 
of V{xi, X2) leads to the following condition: 



^ ^ 25|_,cos2e ' 



The angle is defined as tanfl — v'-^lv'^. The minimization conditions also leads to 

sin2g= (6) 

After B — L breaking, the Zb-l gauge boson acquires a mass ^i'zg^^ — '^g\-iy'^ ■ The high energy 
experimental searches for an extra neutral gauge boson impose lower bounds on this mass. The stringent 



constraint on U{\)b-l obtained from LEP II result, which implies 18| 

My 

> 6 TeV. (7) 

9B-L 

Now we turn to the neutrino sector and show how the observed light-neutrino masses can be obtained 
with 0(1) Dirac neutrino Yukawa coupling. As can be seen from Eq. ([2]), after B — L and EW symmetry 
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breaking, the neutrino Yukawa interaction terms lead to the following mass terms: 



(8) 



where mn = Y^v sin /3 and Mn = Ysv' sin 9. From this Lagrangian, one can easily observe that although the 
lepton number is broken through the spontaneous B — L symmetry breaking, a remnant symmetry: (— 1)^"'"'^ 
is survived, where L is the lepton number and S is the spin. After this global symmetry is broken at much 
lower scale, a mass term for ^2 (and possibly for 5*1 as well) is generated. Therefore, the Lagrangian of 
neutrino masses, in the flavor basis, is given by: 



(9) 



In the basis {i'l, N"^, S2}, the 3x3 neutrino mass matrix of one generation takes the form: 



mo 
m £) Mpf 
A/jv MS2 



(10) 



The mixing matrix(O) for this mass matrix leads to the following light and heavy neutrino masses respectively 
in the limit of << nio, A/jv (where i = 1,2): 



Ml 



Ml 



1 Ml^is, 



2 Ml 



(11) 



where 



O 



y/M%+ml 



3 M%mo^is2 



' V2 4V2 (M^+m2 ) 



mo 



1 1 



3 M%mDp,s^ 
1 MN(M^+4m^)AtS2 



+771% \/2 ^Mj,+ml 4V2 (M^+mJ,)2 V2 ^M^+m|, 4V2 (Af](,+m'^) 



(12) 



Thus, the light neutrino mass can be of order eV, as required by the oscillation data, for a TeV scale Mn, 
provided /isj is sufficiently small, /is^ <^ Mn- In this case, there is no any restriction imposed on the 
value of Dirac mass mu. Therefore, the possibility of testing this type of model in LHC is quite feasible. 
Note that in the limit ~^ which corresponds to the unbroken (—1)^+'^ symmetry, we have massless 
light neutrinos. Therefore, a small non-vanishing 1132 can be considered as a slight breaking of a this global 
symmetry. Hence, according to 't Hooft criteria, the smallness of fis2 is natural. The possibility of generating 
small fis2 radiatively has been discussed in Ref. [3]. 

Finally, it is worth mentioning that the light neutrinos i/i have suppressed mixing (of order mD^iS2/ {Ml + 
iTijy)) with one type of the heavy neutrinos (say fH') and a rather small mixing (of order iriD/Miy) with 
the other type of heavy neutrinos (i^h) by choosing appropriate parameters. The mixing between the heavy 
neutrino and i/'j^ is maximal. In general, the physical neutrino states are given in terms of vl, N'^, and 
^2 as follows: 



1^1 



vh' 



M 



N 



VM_ 
1 



mo 



MNmDfJ.S2 



(M, 



V2 ^Ml + ml 



mo 



V2./MI 



--VL 



+ m|,)3/2 
\/2 



V2 



VMl 
Mn 



V2./MI 
Mn 



S-2 



V2y/Ml 



S2 



-.S2 ^i^L + aiN^ - 02^2 (13) 
~ a{~a2 i^L+ N'- S2) (14) 
a{a2VL+ S2). (15) 
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For niD ^ 100 GeV, Mn 1 TeV and - 1 KeV, one finds that ai ~ 0(10"^°), a2 ~ C(O.l), 
as ^ 0(0.07) and a ^ sin7r/4. In this respect, the gauge eigenstates for neutrinos can be expressed in terms 
of the mass eigenstates as follow: 




( n\ 
\yH' I 



(16) 



III. 



B-L RIGHT-HANDED SNEUTRINO 



In our model, the sneutrino mass matrix of one generation is given by 8 x 8 matrix, which can be decom- 
posed into the following two mass matrices: (i) 6 x 6 mass matrix in basis of {i'L-,i'\,-, N , N'^ , 82-, ■ (ii) 
2x2 mass matrix in basis of {Si,Sl)'^. The ^I's arc decoupled and have now interactions with the SM 
particles. Therefore, one can neglect it and focus on the 6x6 sneutrino mass matrix. In the flavor basis; 
D = {i/L, i>|_, N, N\ S2, Sl)'^ , the sneutrino mass matrix M? is given as 



where 



Mi 



( M^-f - 












\ 






















































Ml.. 




V 








M? - 

S2 S2 





M% . =ml+ cos^ /jy^tye ^ m^Md 



m 



Ml 

^ cos 2/3 ^ 

2 4 



(1 - cot^ 6*), 



M 



Zb-l 



M|t_g^ = m% + \ns2\'^ + MlfMN 



4 

(1 - cot^ 6*), 



(1-cot^ e), 



2 

M^t = fj,*v cos + V sin eVX, MI^.^ = vv' sin 6 sin pY^^V, 
Mj^^ = H'v'cosev^^ +v'smeYX\ Ml^^^= B'^ns^, = ;Us,t;'sin^Y'j. 



(17) 



(18) 

(19) 

(20) 
(21) 
(22) 



In the case of existing the mixing of S2S2 or S2N , the complex scalar DM splits in two real scalar and the 
lighter one is DM. If the mass split is small as well as the momentum of DM, inelastic scattering can be 
considered. 

The mass matrix is diagonalized by unitary matrix F as 



F^M?F = diag(m?m , m?m , m|m , m|m , m|m , m?m ) . 



(23) 



Thus, the mass eigenstates £>™ are defined as Vi = Fiji/™. In general, the lightest sneutrino can be written 
as a linear combination of the sneutrino mass eigenstate, as follows: 



(24) 
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mo [GeV] 

FIG. 1: A contour plot for the mass of B — L right-handed sneutrino in terms of the most relevant parameters: soft 
mass term mo and /i term and for S ~ /? ~ 7r/4 and Yv — 1. The red region represents m^j < 115 GeV, which 
is kinematically excluded because DM mass is lighter than the SM Higgs mass. The green region represents 115 
GeV <msj < 200 GeV. The blue region represents 200 GeV <mp^ < 300 GeV. 



However if one considers large tan/3 and small tan^^ limits, one finds that the lightest sneutrino is mainly 
obtained from the {N, S2) sector. Therefore, it can be expressed as 



t\,N + tIn^ + t\,S2 + T\,Sl (25) 



The mass of this particle, which we call — L right-handed sneutrino" depends on the universal soft 
scalar mass toq and on the parameters and fi'. In Fig. [U we display the mass range oi B — L right-handed 
sneutrino as function of the most relevant parameters mo and fi. If this particle is the lightest SUSY particle, 
then it is stable and can be considered as an interesting candidate for DM. 

Now we consider the relevant interactions of the B — L right-handed sneutrino. From the superpotential 
W in Eq. ^ one gets the following interacting Lagrangian of N in the flavor basis: 



+ Y,%Nl [h,H^2 - iljHt) + Yi.^NfS^ai + h.c, (26) 
Also from the _F-term contributions to the scalar potential one finds the following interaction terms: 

- {Hl^Ht + H^^Hl)Nl{Y''Y'^),flR, ~ iN^Y,DL)ii>lYj N) - ^x* hI\N^Y,vl) 

- {N^YjL){ilY}N) - t,*H^\N^YjL) - iSlY^N){N^YsS2) - fi*xUN^YsS2) 

- \xi\'n!{YsY^)N, - fi*s,XiNlYs,jSl^ + h.c. (27) 

Next the interactions of N with the gauge fields lead to the following Lagrangian: 

^?nt = 9l-LZl-LN!N^ - Wb-lZ^s_^ (^.N^ ~ N^d^N,) - iV2gB-LN,N^ZB-L + h.c. (28) 
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Finally, the _D-term implies that 



X 



(29) 

In our analysis for B — L right-handed sneutrino {vi) DM, the following assumptions are considered for 
simplification: (i) Each sector consists of one generation only, (ii) nisM < ^dm < "m-susY, fnfi±, m^-^^, 
Mzb^lj "^Si 2 1 where msu is standard model particles including of the lightest neutral Higgs boson (/iq), 
msuSY is the supersymmetric particles, and itidm is the mass of the lightest SUSY sneutrinos. {Hi) In 
chargino sector, mass is approximately given by /x, in the limit of /x, M2 » raw, where mw is the 
SM charged weak gauge boson mass. Under these assumptions, the relevant interacting Lagrangian of vi 
consists of 

Note that the D-term is kinematically irrelevant, while the F, G interactions take the following simplified 
form: 



h.c, (31) 



£L = -iYjY,)T,^T3^\hyi>l^i, (32) 

^?nt = -■i9B-LZ'^_LT4iT3i (l>idf,l)l - i>l9^I>i) , (33) 

where is considered to be included in the mixing of Higgsino and Chargino. Hereafter we use it as Y^. 
Eq. ([3T|) might be relevant to the indirect detection, which will be discussed in the section IVll Eq. ([32]) is 
more relevant to the WMAP experiment, which will be discussed in the next section. Eq. (|33p is applied to 
analyze the direct detection as CDMS II/XENON, which will be also discussed in the section fVl 



IV. RELIC ABUNDANCE OF B - L RIGHT-HANDED SNEUTRINO 

In this section, we compute the relic abundance oi B — L right-handed sneutrino DM. We consider the 
standard computation of the cosmological abundance, where I'l is assumed to be in thermal equilibrium with 
the SM particles in the early universe and decoupled when it was non-relativistic. Therefore, the i>i density 
can be obtained by solving the Boltzmann equation: 

^ + 3Hn,, = < atrv > [(n.J^ - (n?^')'], (34) 

where is i>i number density with p/^-^ = mj^-^rij^-^. One usually defines flj^-^ = pc^/ Pc, where Pc is the critical 
mass density. In addition, < cr?""7j > is the thermal averaged of the total cross section for vi annihilation 
into SM lighter particles times the DM relative velocity v. For non-relativistic f^i, the thermal averaged 



annihilation cross section, < cr-""?; >, can be approximated as follows |22| : 

<<"%>~a5, (35) 
where a^i bo^ are the coefficients coming from s-wave and p-wave of vivi annihilation, respectively. 



8 



1^1 



1^1 



9+ 



FIG. 2: Possible annihilation channels of v\. The second diagram gives a sub-dominant contribution, however it may 
be relevant for indirect detection processes. 



From Eas. pm55)) . one finds that the dominant annihilation channels of vi are given in Fig. [2] It turns 
out the annihilation oi B — L right-handed sneutrino into SM-like neutral Higgs, through the four point 
interaction vertex, gives the dominant contribution. The tree level annihilation channel DiVi — > £^£~ is 
suppressed by the mass of the chargino exchanged particle. This channel may be relevant for the indirect 
detection processes which will be discussed later. Our computation for the annihilation cross section leads 
to the following ap^ bj^-^ : 

= ^/^l>'-^3l^4l|^ (36) 

1^1 

where 

= m,, ' '^'^ ■^'^ 2(1 - zl) 



^, p'^^^l-zl:,l = ,^^^^^, (38) 

2 



here we define that a is a final-state particle. From Eg. ([Ml) one finds that the relic abundance ^c,^h^ is given 



by m 



where 



0.0955 mpi TOpi(ajii 4- e^iii/xF) .,„^ 
= m -J- ■ . (40) 

Here m-p\ is the Planck mass (1.22 x 10^^ GeV) and gt{Tp) enumerates the degrees of freedom of relativistic 

1/2 

particles at the freeze out temperature Tp, which can be fixed as gj (TV) = 10. From the above expressions, 
one notes that vi relic abundance depends only on the right-sneutrino mass and the annihilation cross section 
coefficients ap^ and ba-^ . 

Given ap^ bci^ we can determine the freeze out temperature Tp, below which the B — L right-handed 
sneutrino annihilation rate is smaller than the expansion rate of the universe and then computing the relic 
density fiij^/i^. In our numerical computation we assume a universal soft mass (toq) and fix the other 
parameter as follows: 

niho = 115 GeV, Afz^-i, = 6 TeV, Mn = 1 TeV, ^isi,2 = 1 KeV, v = 175 GeV, 

m£, = 100 GeV, 2 = i'' = M*'' = 1 TeV, 0.1, F'^ = 0.05, i;-^ = = 0.1 TeV, (41) 
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FIG. 3: Relic density flcj/i^ as function of the B — L right-handed sneutrino mass m^j for three values of Dirac 
neutrino Yukawa coupling. The region between two lines is allowed by the experiment of WMAP P]. 

In addition we analyze the relic density of vi in the following regions of the parameter fi, and 9: 

< 6* , ^ < TT, 50 GeV < , toq < 300 GeV. (42) 

In Fig. [3l we present the values of relic abundance oi B — L right-handed sneutrino floih^ as a function 
of for the following values of Dirac neutrino Yukawa coupling: YJy = 0.5,0.7, 1.0. We require B — L 
right-handed sneutrino relic density to be 0.09 < Vla^h^ < 0.12 in order to be consistent with WMAP results 
at 3(7 2J|. As can be seen from this figure, smaller values of Dirac neutrino Yukawa coupling are favored 
and lead to more allowed points that satisfy the WMAP observational limits of DM relic density. Also we 
find that the smaller Y^, is considered, the smaller DM mass one obtains. 



DIRECT DETECTION 



In this section we discuss the possibility to detect our B — L right-handed sneutrino in direct detection 



experiments such as CDMS (II) ;2] and XENON 100 experiment 
scattering cross section between DM j>i and nuclei N is given by 

t2 



22 



The reduced mass 



is defined as 



lOTT 



1 



1 

M 



The general form of the elastic 



(43) 



(44) 



9B-L 



Zb- 



L 



q 9B-L q 

FIG. 4: Our dominant diagram for the direct detection 
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Mz /gB.L>6.0TeV 
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XENON100 




100 1000 

[GeV] 

FIG. 5: Elastic scattering cross section between a nuclei and the B — L right-handed sneutrino DM as function of 
the DM mass for Qb-l ~ 1- 



where M is the nuclei mass. The coefficient is given by 



9N 



{A^Z)b„ + Zbp, 



2b„ 



2b, 



(45) 



Here A and Z are the mass number and the atomic number, respectively. The effective Lagrangian parameters 
bu and bd are defined as 



£eff = bgXt'qjt'q, q^{u, d). 



(46) 



Here is a general form of the vector current. In case of fermionic DM is given by ~ fl^f- While 
for bosonic DM, it is defined as X'^ ~ ib^d'^b — ibd^b'^ . 

In our B ~ L case, the elastic scattering cross section of the right-handed sneutrino with a given nuclei 
has a spin-independent contribution arising from Zb-l gauge boson exchange diagrams, as can been seen 
in Fig. m The interactions between the right-handed sneutrino vi and Zb-l boson in B — L model are 

1 



C D -igB-LZ'^_j^TiiT3i [I'ldfj^i'l - v\df^vi 
Writing down the effective interaction, one obtains 

|r4ir3i|2g|_i 



1 



iCcS 3 -- 



Assuming |r4ir3ip ~ 1 for simplicity, one finds 



bp = b„ 



Therefore, b^ is given by 



ON 



lA 



■ 9%-L 

'Ml ■ 



9%-L 

Ml 



Thus, the elastic scattering cross section oi B — L right-handed sneutrino is given by 

Zb-l _ "^r \a\2_9b~l 



(47) 



(48) 



(49) 



(50) 



(51) 
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In Fig. Owe depict the relation between the cross section and B — L right-sneutrino mass. As can been 
seen from this figure, the following upper bound on cr?_I^S^ can be obtained: 

ct~,''Sn < 6-2 X lO-''^ cm^. (52) 

It is also remarkable that the elastic cross section is quite insensitive to the B — L right-handed sneutrino 
mass mpj. However, one observes that a light sneutrino ^ 0(100) GeV is more favored by direct detection 
experimental results. The current limits from CDMS II and XENON experiments indicate to a lower-bound 
of order 3.7 x lO"**"* cm^. This suggests that our B~ L right-handed sneutrino DM is expected to be detected 
in the direct detection experiments in near future. 

Before concluding this section, it is worth mentioning that the XENON 100 [2^ experiment has recently 
presented new limits on the WIMP-nucleon cross section for inelastic DM. These limits are due to a DM run 
with 100.9 live-days of data, taken from January to June 2010. It was shown that (J^-n < 10^**^ cm^ can 
be extracted for rriy^ > 100 GeV. Where x is a generic dark matter. The bound rules out the explanation of 
controversial DAMA/LIBRA modulation results, as being due to inelastic DM. The bound obtained on o-^^n 
from the inelastic DM analysis of XENON 100 should be considered carefully, since a minimum velocity for 
DM to scatter in a detector is introduced, hence a large amount of fiducial is needed. Nevertheless, as can 
be seen from Eq. (52), the cross section of our B — L sneutrino is well below the XENON 100 bound. 



VI. INDIRECT DETECTION 
A. PAMELA and Fermi-LAT experiments 



34 



As advocated in the introduction, the indirect searches for DM by the Space Observatory PAMELA 
and Fermi-LAT 35], indicate that the DM may contribute to the positron flux by direct annihilation in 
This is one of the main feature of our DM candidate B — L right-handed sneutrino, therefore, it 
important to investigate the possibility that it accounts for these results. PAMELA collaboration reported 
excess flux between 8 and 80 GeV, with no excess in the corresponding anti-proton flux. Also ATIC and 
Fermi-LAT balloon experiments have shown excess electron and positron flux at energies around 10 — 1000 
GeV. While there are plausible astrophysical explanations for these excesses, such as local pulsars and 
supernovas remnants, they could also result from DM annihilation. Note that in order to explain the Fermi- 
LAT experiment by DM annihilation, the DM mass must be of order 0{TeV). As shown in the previous 
section, such heavy DM mass is not favored by of direct detections. Therefore, in our analysis, we assume 
that the Fermi-LAT experiment may be saturated by considering an astrophysical background and we will 
focus on PAMELA measurement. 

It is known that if the DM annihilation is to explain the observed anomalous flux, a large annihilation cross 
section, (av) ~ 10"^'' cm^s""'^, is required to fit the excess flux, which is incompatible with straightforward 
estimates of the relic DM abundance in conventional cosmological models. Otherwise, a huge, unexplained, 
boost factor must be introduced. In our SUSY B ^ L model, the B ~ L right-handed sneutrino annihilates 
into t^i~ channels, as shown in the second diagram of Fig. [2l However, as discussed in the previous section 
that these channels give sub-dominate contribution to the annihilation process. Therefore, the corresponding 
annihilation cross section is < 10^^^ cm'^s^^. In this case, one requires a huge boost factor O(10^~^) at 
least in order to account for PAMELA results. In general, it is known that there are two mechanisms to 
enhance the cross section and may justify this large boost factor. The first is Breit-Wigner mechanism [36 [ 



and the second is the Sommer-feld [33] enhancement. Breit-Wigner mechanism can not be implemented in 
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our model, since there is no any diagrams with s-channeL On the other hand, the Sommer-feld enhancement 
requires higher DM mass compared to the mediated particles in order to obtain enough large boost factor 
> 0(10^"^). This assumption we have already avoided in order to not spoil the direct detection. As a 
result, it is difficult for our B — L right-handed neutrino to explain the controversial results of PAMELA 
experiment. 



B. Muon flux measurement from Super-Kamiokande 



The high energy neutrinos induced by DM annihilations in the earth, the sun, and the galactic center is 
an important signal for the indirect detection of DM. Such energetic neutrinos induce upward through-going 
muons from charged current interactions provide the most effective sig natures in Super-Kamiokande. The 
neutrino-induced muon flux is evaluated from the neutrino flux 



40 



41| as 



F 



(ann) 



5.9 X 10 



cm ^s ^ X 



{av)F 



(J2)oA» 
10 



(53) 



where we fix the typical values (J2)oAr2 ^ 10 for V'max = 5° in case of the Navarro-Frenk- White (NFW) 
halo density profile [42j, and F collectively denotes the primary annihilation mode (e.g., r+r^, etc.). Notice 
that model dependence comes from {av)p, which will be shown later. Sp is defined as 



dN 



dE 



P \ — 



dE, 



(54) 



where Ei^ = W£;^ , and iJmin is the threshold energy above which the muons can be detected. -Pi/ii/^ denotes 
the probability that the i^i at the production is observed as i/^ at the Earth due to the effect of neutrino 
oscillation. Regardless to the complicated expression of Sp, it is found as a fixed value depending on each 
of the F particle, which is, e.g., 0.2 for fj, pair, 0.14 — 0.18 for t pair, 0.78 for Vt- pair, etc js^. The limits 
from Super-Kamiokande are given in the ref. 43 1. 



Ft 



1 



i+ 



FIG. 6: Fenymen diagram of uii'i — !■ £ . Note that, in general, there is an another contribution of vi^ui — >■ hP 
V, V process. However, due to the masslless limit of v, the cross section of this channel vanishes. 



In SUSY B — L model, the relevant Lagrangian in given by 

£Zt = Y,Tiiu\ W-Plul - h+^PL^l + h.c. 



(55) 
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One can show that in the hmit w — > 0, the thernLahy averaged cross section is given by 



{av)F\v^o 



y 



|i;r3ir4i|^4(2 - zl) X + rlj^_ + rpj^^Tp -^^ + rlj^^ 



E 

F=T 



167rm| 



|l^r3ir4i|S|(2-z|,) X r 



(56) 



Note that the term proportional to ?'^^o vanishes in the massless hmit of v. The parameters zp and P'p are 
as defined in Eq. ((38|) . The other parameters are defined as foUows: 



Wo 



rpo 



{l^zl + wl)-\ 



(57) 



Here F refers to the final-state particle and a denotes the mediated-particle. In our analysis for the muon 
fiux, we consider the same set of inputs that we have used in the previous sections that leads to relic 
abundance within the WMAP limits: 0.09 < Uh^ < 0.12. Moreover, since we are considering the effect one 
generation, we assume r final state only. Using Sf that corresponds to r final state, one finds (av) 



F 



(ann) 



5.9 X 10^15 cm-^s-i X J2 Sf 

F = T± 



{(JV)F 



6.9 X 10"^ cm^^s^i X 



F = T± 



3.56 X 10-7 



(J2)oA» 
10 



167rm?^ 



i^r3ir4ir4(2-z:^) x 



F.h+ 



(58) 



As we fixed the cone-half angle from the galactic center: ^/imax = 5° which is maximum in the case of NFW 
profile, one finds that the Super-Kamiokande limit should be less than 5 x 10-^''^cm-^sec~^ (See figures of 
Ref. 



38|). 



In Fig. [71 we plot the muon flux induced from the B ~ L annihilation in the galactic center. As can 
be seen from this figure, the result of the induced muon flux is < lO-^*' cm~^sec~"'^, which is few order of 
magnitudes smaller than the upper bound of Super-Kamiokande. This result is expected since the thermally 
averaged cross cross is much smaller than the typical value required by indirect detections. 



r-.(ann) r —2 — li 

F\ ' cm sec 



10 



10 



10 



10 



10 



10 



-10 
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mc, [GeV] 



250 300 



FIG. 7: Muon flux induced from the B ~ L annihilation in the galactic center. The horizontal line refers to the 
Super-Kamiokande upper bound which is given by ~ 5 x 10~^^cm~^sec~^. 
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VII. CONCLUSIONS 

In this paper we considered the supersymmetric B — L model with inverse seesaw mechanism. We demon- 
strated that the hghtest right-handed sneutrino in this model can be stable and a viable candidate for cold 
dark matter. We studied the relic abundance of the B — L right-handed sneutrino and showed that the 
WMAP result, ilh"^ ~ 0.11, can be satisfied in a wide range of the parameter space. We emphasized that 
the dominate annihilation channel oi B — L right-handed sneutrino is given by vivi — ?> hPh^ , where hP is the 
SM-like Higgs. We also studied the direct detection rate oi B — L right-handed sneutrino. We found that 
its elastic cross section is consistent with the upper bounds of current experiments, such as CDMS II and 
XENON. Our result oi B — L right-handed sneutrino direct detection is promising and indicates that it can 
be detectable in near future experiment. 

In addition, we have analyzed the indirect detection rate oi B — L right-handed sneutrino. In particular, 
we focused on the observation of the Space Observatory PAMELA for positron flux and also on the neutrino- 
induced upward through-going muons in the Super-Kamiokande detector. We showed that although the B—L 
right-handed sneutrino annihilates at tree level into leptons, the corresponding cross section is much smaller 
than the required one for accommodating PAMELA results. Also the neutrino flux induced by the B — L 
right-handed sneutrino annihilations in the galactic center is much smaller than the Super-Kamiokande's 
limits. 

Finally it is worth mentioning that our B — L right-handed sneutrino may be produced at the Large 
Hadron Collider (LHC) through the channel is qq — > Zb-l — > i^li'i- However, the amplitude of this channel 
vanishes identically due to the fact that left and right quarks or leptons has the same B — L quantum 
numbers. However, slepton/left-handed sneutrino may decay to right-handed sneutrino, which escapes the 
detector and gives the missing energy signal similar to other examples of cold DM j46i] . 
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